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1.  SUMMARY 


An  optimal  design  procedure  of  the  composite  flywheel  rotor  is  proposed  with  a 
consideration  of  the  residual  stresses  caused  during  the  manufacturing  process.  The 
composite  materials  are  inherently  anisotropic:  the  mechanical  and  thermal  properties 
differ  between  the  longitudinal  and  the  transverse  directions.  This  anisotropic  property 
makes  the  composite  materials  superior  to  other  conventional  metal  in  the  design  of  the 
flywheel  rotor.  However,  at  the  same  time,  the  anisotropic  properties  cause  the  residual 
stresses  to  develop  in  manufacturing  process.  The  induced  residual  stresses,  especially  in 
the  weak  radial  direction,  deteriorate  the  overall  performance  of  the  composite  flywheel 
rotor.  The  residual  stresses  in  the  thick  wound  composite  structure  can  be  as  high  as  the 
strengths  of  the  rotor.  Such  residual  stresses  add  to  the  mechanical  stresses  developed 
during  the  rotational  speed,  thus  limiting  its  full  capability.  Therefore,  in  addition  to  the 
mechanical  stresses,  the  residual  stresses  should  be  controlled  in  the  manufacturing 
process  and  properly  taken  care  of  in  the  design  of  the  rotor.  The  thorough  simulation  of 
the  manufacturing  process  and  the  residual  stresses  is  extremely  difficult  since  many 
physical  phenomena  are  involved.  In  this  work,  first  of  all,  a  simple  model  has  been 
proposed  which  quantitatively  characterize  the  process-induced  residual  stresses  using  a 
effective  temperature  change.  The  cylinders  of  three  different  sizes  have  been 
manufactured  and  the  residual  strains  in  the  radial  and  hoop  directions  are  then  measured 
by  split-ring  method.  The  measured  residual  stains  have  then  been  predicted  and 
compared  with  those  by  the  model.  The  effects  of  the  residual  stresses  on  the  hybrid  rotor 
were  then  analyzed  by  comparing  the  two-dimensional  solution  to  the  three-dimensional 
finite  element  program.  The  residual  stresses  in  the  rotor  with  the  angle  variation  along 
the  radial  direction  were  also  analyzed.  Finally,  the  stress  analysis  module  was  linked  to 
an  optimization  program  in  which  the  design  variables,  e.g.,  the  thickness  of  each  hybrid 
ring  can  be  optimized  for  ^the  various  size  and  material  sequences  of  the  composite 
flywheel  rotor  with  a  consideration  of  both  the  residual  stresses  and  the  rotational 
centrifugal  stresses. 
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2.  INTRODUCTION 


An  advanced  composite  flywheel  has  recently  attracted  the  considerable  attention 
of  many  investigators  for  many  energy  storage  applications  including  electric  utilities, 
hybrid  or  electric  vehicles  and  spacecraft  [1-5].  The  composite  flywheel  rotor  has 
characteristics  of  distinctively  high  energy  density,  long  life  and  lightweight.  The 
essential  components  of  a  flywheel  energy  storage  system  are  a  composite  flywheel  rotor, 
motor/generator,  power  transforming  electronics,  magnetic  bearing,  and  housing.  Recent 
efforts  in  development  of  the  flywheel  have  been  devoted  to  the  experimental  and 
theoretical  research  on  increasing  the  performance  and  insuring  the  safety  of  the  flywheel. 
It  is  well  known  that  the  rotor  strength  can  be  significantly  improved  by  properly 
adjusting  the  ply  stiffness  along  the  radial  direction.  Especially,  along  with  technical 
innovations  in  high  strength  composite  materials,  many  investigations  have  been  devoted 
to  increasing  the  total  stored  energy  (TSE)  of  the  flywheel  rotor  [6-10].  The  approaches 
for  calculating  the  stress  distribution  of  a  hybrid  composite  rotor  subjected  to  centrifugal 
forces  have  been,  cited  in  literature  [1-10]. 

During  the  curing  process  of  the  filament  wound  structure,  the  residual  stresses  are 

developed  mostly  due  to  the  chemical  shrinkage  and  the  thermal  expansion.  Those 

residual  stresses  are  very  complicatedly  affected  by  the  degree  of  anisotropy,  the 

dimensions  of  the  rotor  and  the  curing  temperature  cycle  [12-16].  The  radial  tensile 

stresses  developed  during  the  curing  process  generally  deteriorate  the  energy  storage 

capability  since  the  stresses  are  additive  to  the  radial  tensile  stresses  which  are  also 

developed  during  the  rotation  of  the  rotor.  Therefore,  the  curing  stresses  have  to  be 

reduced  as  much  as  possible  to  increase  the  energy  storage  capability.  However,  the 

certain  level  of  the  curing  stresses  is  indispensable  and  cannot  be  completely  avoided. 

Therefore,  they  have  to  be  considered  in  the  design  of  the  composite  flywheel  rotor. 

However,  methods  for  evaluating  the  effects  of  material  properties  and  curing  residual 

\ 

stress  on  the  energy  storage  ^capability  have  not  been  pertinently  addressed.  Moreover,  no 
investigation  thus  far  has  addressed  how  to  optimally  design  the  hybrid  composite  rotor 
considering  both  the  curing  residual  stresses  and  the  centrifugal  stresses. 

In  this  work,  the  residual  stresses  developed  during  the  curing  process  are  modeled 
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as  an  effective  thermal  loading,  which  greatly  simplify  the  calculation  of  the  residual 
stresses;  in  this  way,  only  a  temperature  increase  need  to  be  determined.  The  two- 
dimensional  stress  analyses  for  the  flywheel  rotor  based  on  the  assumptions  of  plane 
stress  (PSS),  plane  strain  (PSN),  and  generalized  plane  strain  (GPS)  have  been  developed 
which  consider  the  material  hybridization  and  the  angle  variation  together  with  the 
thermal  loading.  The  two-dimensional  analyses  are  then  evaluated  using  three 
dimensional  finite  element  methods.  Overall  it  is  found  that  PSS  is  the  best  analysis 
method  for  the  flywheel  rotor.  The  equivalent  thermal  loading  is  also  estimated  using  the 
well-known  cure  simulation  program  WINDTHICK  [14-16],  and  it  is  found  that  the 
method  is  very  accurate  enough  to  predict  the  residual  stresses  and  strains  developed 
during  the  whole  curing  process. 

An  optimal  design  program  (FLYOPT)  has  been  finally  developed  using  the 
Microsoft  Excel  interface.  In  the  program,  the  rotor  is  divided  into  as  many  rings  as  users 
define  and  each  thickness  and  winding  angle  of  the  rings  can  be  optimized  to  maximize 
the  performance  of  the  flywheel  rotor,  especially  the  total  stored  energy. 


3.  2-DIMENSIONAL  STRESS  &  STRENGTH  ANALYSIS 


For  a  numerical  calculation,  the  composite  flywheel  rotor  of  varying  material 
properties  is  divided  into  many  rings  and  each  ring  is  assumed  to  have  a  constant  material 
property  of  an  angle-ply  laminate  [(±<p)m]  as  shown  in  Figure  1.  Based  on  the 
assumption  of  a  generalized  plane  strain  state,  all  the  rings  are  free  to  expand  or  contract 
in  the  axial  direction  but  have  the  same  axial  strain.  The  flywheel  rotor  is  also  assumed 
cylindrically  orthotropic  and  subject  to  axisymmetric  centrifugal  forces  due  to  rotation. 
With  the  axisymmetric  conditions,  the  circumferential  displacements  are  thus  neglected. 
All  the  displacements  and  stresses  are  thus  independent  of  the  axial  and  circumferential 
directions.  The  radial  equilibrium  equation  for  each  ring  and  the  axial  force  equilibrium 
equation  for  the  rotor  are  concurrently  solved  and  expressed  in  a  stiffness  matrix  with  the 
radial  displacements  and  the  uniform  axial  strain.  Once  the  displacements  and  the  axial 
strain  are  obtained,  the  stresses  and  strength  ratios  are  calculated  for  all  the  rings. 

3.1  Stress  Analysis 


The  stress  distribution  in  each  ring  is  governed  by  the  radial  equilibrium  equation, 
which  is  written  in  cylindrical  coordinates  as  [24] 


d(J  (7  (7 a  2  r\ 

— -+  — - -  +  pro)2  =0 

dr  r 


(i) 


where  ar  and  Gq  are  the  radial  and  circumferential  stresses,  respectively,  p  denotes  a 
density,  and  co  the  rotational  angular  velocity.  The  stress-strain  relation  in  the  cylindrical 
coordinate  system  as  shown  in  Figure  1  is  written  as 


(7-Q(£-  aAT  )-Qe-  fiA  T 


(2) 


where  crand  fare  the  stress  and  strain  vectors,  respectively,  Q  and  a  denotes 
matrix  and  thermal  expansion  coefficient  vector  in  cylindrical  coordinates.  Eq. 
rewritten  in  the  cylindrical  coordinate  system  as 
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In  Eq.  (3),  the  axial  stress  is  zero  in  a  plane  stress  (PSS)  assumption  and  the  axial  strain  is 
zero  in  a  plane  strain  (PSN)  assumption.  More  generally,  the  axial  strain  is  assumed  to  be 
a  constant  in  the  generalized  plane  strain  (GPS)  assumption.  Together  with  these 
assumptions,  the  circumferential  and  radial  strains  are  expressed  by  the  r-direction 
displacement  ur : 


£d 


u  du 

r  or 


=  0,  and  \ 


<JZ  =  0  PSS 
ez  =  0  PSN 
£z=£a  GPS 


(4) 


Substituting  Eqs.  (3)  and  (4)  into  Eq.  (1)  yields  the  radial  displacement  ur: 

ur  =  -pa1  <p0r}  +  c,  q\rK  +  C2  <p2r~K  +  <p3eor  +  <pTIA  Tr  (5) 

C]  and  C2  in  Eq.  (5)  are  unknown  constants  to  be  determined  from  the  boundary 
conditions,  and  K*and  (pi  are  defined  in  terms  of  the  material  properties: 

1  1 


K=  ^ 


Qr/  r"  (9->r)Q33,<P'  Q„  +  kQ„ 

- 1 - ,  <p  =  On  Z&l  and  (p  =AzA_ 

0.3  "*033  033-01.  033 ~Qu 


<Pl  = 


Substitution  of  Eq.  (5)  into  Eqs.  (3)  and  (4)  yields 

£r  =  - 3pco2<por 2  +  kCx(P/k~]  -  KC2(p2r~K~l  +  <p2£0  +  cpnAT 
Gr  ~  - pco2(pAr 2  +  C,rA_1  +  C2r~K~]  +  (ps£o  +  <pT2AT 

where 


(6) 

(7) 

(8) 
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In  order  to  derive  the  stiffness  matrix,  the  radial  displacements  at  the  inner  and 
outer  surfaces  are  written  in  a  displacement  vector  using  Eq.  (4), 


U  =  pCD2Ua)Jr£0U£+&TuT+G&C 


(10) 
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In  Eq.  (11),  r,  denotes  the  inner  radius  of  a  ring,  and  r0  the  outer  radius  of  a  ring.  The 
normal  radial  stresses  at  those  two  radii  can  be  also  written  in  a  force  vector  using  Eq.  (8), 
ft  =  pco2fa+£0f£+ATfT+fGC  (12) 
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Eliminating  the  constants  C  in  Eqs.  (10)  and  (12)  derives  the  relation  of  stresses  and 
displacements  at  the  inner  and  outer  surface  of  each  ring.  The  derived  equation  can  be 
expressed  as 


ku  =  fb  +  pco 2fa  +  £jc  +  ATfT 


The  ring  stiffness  matrix  k  is 


k  =  l‘G0'G =  — 

<r2 


CiQ\3  -  xQr- 
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(14) 


(15) 


where 


c,=rr+r,f2=rr-rand^i  06> 

ro 

Notice  that,  in  defining  the  force  term  ft,  of  Eq.  (14),  the  stress  component  was  multiplied 
by  its  radius,  and  the  ring  stiffness  matrix  k  becomes  symmetric.  In  Eq.  (14),  fa*  ft  and 
fT  are  related  to  the  centrifugal  force  and  defined  as 


4  = 
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The  constant  axial  strain  £0  in  Eq.  (14)  can  be  computed  from  force  balance  in  the 
axial  direction.  Since  there  is  no  external  axial  force  applied  on  the  top  and  bottom 
surfaces,  the  summation  of  the  axial  forces  over  all  the  rings  should  be  zero: 

Y,27iti'<7^rdr=Q  og) 

7=1 


Using  the  expressions  of  oz  in  Eq.  (3),  and  ur  and  £r  in  Eqs.  (5)  and  (7),  Eq.  (18)  is 
expressed  as 


In  Eq.  (19),  j  denotes  a  ring  number  and  ^  ofa  y-th  ring  is  defined  as  follows  (the 
superscript  j  is  omitted  for  clarity): 


where 
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In  Eq.  (19)  and  Eq.  (21),  the  following  symbols  are  defined: 
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(22) 


Eqs.  (1)-(19)  apply  to  any  arbitrary  ring.  At  each  interface  between  two  adjacent 
rings  (e.g.,y‘  and  y+ 1),  the  continuous  radial  traction  and  continuous  radial  displacement 
should  be  satisfied  as  illustrated  in  Figure  2: 

<7<y+1)  =  cr</)  and  «<y+1)  =  (23) 

Combining  all  the  ring  stiffness  matrix  in  Eq.  (14)  by  using  the  continuity  equations 
yields  the  global  stiffness  matrix  equation 
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(24) 
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The  global  displacement  vector  U  consists  of  N+\  displacements  at  all  the  radii  of  the 
rings.  The  global  force  vectors  F/,,F  a ,  F'^  F^and  FV  in  Eq.  (24)  are  as  follows; 
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7=1  7=1  V  7 
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Notice  that  the  global  stiffness  matrix  is  symmetric,  and  F*  has  only  the  boundary 
pressures  given  at  the  inner  and  outer  surfaces  since  the  internal  stresses  were  canceled 
out  due  to  the  stress  continuity  condition.  F<y  is  an  internal  inertial  force  vector  due  to  the 
rotation. 

Now  we  can  solve  the  system  of  7V+2  equations  and  obtain  the  radial  displacements 
at  each  interface  and  constant  axial  strain  X  The  strains  in  cylindrical  coordinates  are 
then  calculated  using  the  relationship  in  Eq.  (4). 


3.2  Strength  Analysis 


The  on-axis  strains  e  in  the  material  symmetric  axes  are  obtained  by  transforming 
the  off-axis  strains  £  in  the  cylindrical  coordinates[25]  as  shown  in  Figure  2: 

e  =  Ts  (27) 

where  T  is  a  transformation  matrix  defined  in  terms  of  a  ply  angle  (j) ; 
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The  on-axis  stress  a  is  then  calculated  using  the  stress-strain  relationship: 

a-Q{e-aAT)  (29) 

The  on-axis  stiffness  matrix  Q  and  the  coefficients  of  thermal  expansion  (CTE)  a  in 
Eq.  (29)  are  given  by  the  material  properties [2 5]: 
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The  definitions  of  Moduli,  Poisson’s  Ratios,  and  CTE  in  Eq.  (30)  are  explained  in  Table  1. 

Once  the  on-axis  stresses  in  the y-th  ring  are  calculated,  a  failure  of  the  material  can 
be  accessed  by  using  a  failure  criterion.  In  this  study,  3-dimensional  Tsai-Wu  quadratic 
failure  criterion[19]  is  used: 

aTFa+FaR-R2=0  (31) 

Solving  Eq.  (31)  yields  a  strength  ratio  R  indicating  material  failure  with  a  value  greater 
than  1 .  When  the  strength  ratio  R  <  1  ,  the  applied  stresses  can  increase  by  a  factor  of  its 

inverse  before  material  failure  occurs.  The  strength  parameters  F  and  F  are  defined 
in  terms  of  the  material  strengths: 


jF  = 


1 

111 

111 

XX' 

2\XX'YY’ 

2  V  XX'  YY' 

i  , 

1  1 

1 

1 

2  V 

XX' YY' 

YY' 

2YY' 

I  , 

1  1 

1 

1 

2  i 

XX'  YY' 

2  YY' 

YY' 

0 

0 

0 

s2J 


(32) 


F  = 


1 _ 1_ 


I__L 

Y  Y' 


1-1  0 

Y  Y 


(33) 


where  X  and  X’  are  the  tensile  and  compressive  strengths  in  the  fiber  direction,  Y  and  Y1 
those  in  the  matrix  direction,  and  S  the  shear  strength. 
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4.  COMPARISON  OF  2-D  &  3-D  SOLUTIONS 


Three  types  of  two-dimensional  analyses  based  on  the  assumptions  in  the  axial 
deformation  and  stresses,  i.e.,  plane  stress  (PSS),  plane  strain  (PSN),  and  generalized 
plane  strain  (GPS)  are  described  in  the  previous  chapter.  The  accuracy  of  the  two- 
dimensional  analyses  of  the  flywheel  rotor  with  the  residual  stresses  depend  upon  the 
length  and  material  properties  of  the  rotor  and  the  boundary  condition  at  the  axial  ends  of 
the  rotor.  However,  either  one  of  them  cannot  perfectly  match  with  three-dimensional 
analysis.  In  this  chapter,  the  two-dimensional  analyses  will  be  compared  with  the  three- 
dimensional  analysis  using  the  commercial  finite  element  analysis  program  ANSYS.  The 
full  three-dimensional  solutions  with  a  consideration  of  the  rotor  length,  the  material 
hybridization,  and  the  helical  winding  are  obtained  and  compared  with  three  different 
two-dimensional  solutions,  i.e.,  PSS,  PSN,  and  GPS.  The  symmetry  of  the  rotor  enables 
to  analyze  the  one-eighth  portion  of  the  rotor  in  the  three-dimensional  finite  element 
analysis  of  the  rotor  as  shown  in  Figure  2.  Eight-node  hexagonal  element  with  each  node 
having  three  degree-of-freedoms  is  used  for  the  analysis. 

Total  of  6  cases  are  considered,  all  with  the  same  ID  as  75  mm  and  OD  as  221  mm. 
First  3  cases  are  for  the  hybrid  rotor  of  Glass/Epoxy  and  T300/5208,  and  the  rest  3  cases 
are  for  the  helically  wound  rotor  of  T300/5208.  The  specifications  of  the  cases  are 
summarized  and  listed  in  Table  1.  Their  material  properties  are  shown  in  Table  2.  The 
case  1  is  a  reference  case  to  compare  with  other  cases  of  hybrid  rotors.  In  case  2,  the 
length  is  changed;  in  case  3,  the  rotational  forces  are  applied  instead  of  the  thermal 
loading.  In  case  4,  the  winding  with  the  angle  of  45  degrees  for  the  half  of  the  total  radial 
thickness  and  then  hoop  winding  for  the  rest  of  half  are  considered.  In  case  5,  the  length 
is  changed,  and  in  case  6  the  helical  winding  angles  are  linearly  changed  from  45  degrees 
at  the  inner  surface  to  0  degree  at  the  outer  surface. 

The  radial,  hoop,  and  axial  stresses  and  strains  for  the  6  cases  are  calculated  by  both 
three-dimensional  finite  element’ and  the  two-dimensional  methods,  and  they  are  shown 
in  Figures  3~8.  In  the  three-dimensional  analyses,  the  solutions  are  obtained  at  the  middle 
and  top  surfaces  of  the  rotor.  The  stresses  and  strains  at  the  top  surface  in  the  longer  rotor 
slightly  but  negligibly  differ  from  those  of  shorter  rotor,  as  shown  Figures  3  and  4.  The 
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PSS  and  GPS  solutions  match  well  with  the  3-D  solutions  at  the  middle  surface,  which 
can  be  considered  to  be  representative  of  the  rotor  status  since  they  show  a  little  higher 
values  than  at  the  top  surface.  On  the  other  hand,  the  PSN  solution  is  far  higher  than  other 
solutions.  Notice  that  the  axial  stains  can  be  well  predicted  using  the  PSS  solution.  In  the 
case  3,  both  GPS  and  PSS  solutions  also  well  predict  the  3-D  solution.  The  axial  strain  in 
case  3  can  be  best  predicted  by  the  PSS.  In  cases  4  and  5,  the  stresses  at  the  top  are  quite 
higher  than  those  at  the  middle.  We  can  notice  that  the  abrupt  change  of  winding  angles 
from  45  to  0  degree  causes  the  high  radial  stress  variation  around  the  interface,  which 
cannot  be  well  predicted  by  any  of  2-dimensional  analyses.  The  radial  stress  variation 
around  the  interface  at  the  top  surface  is  higher  in  the  case  5  in  which  the  rotor  is  longer 
than  the  case  4.  However,  except  the  region  near  the  top  surface  the  GPS  solution  is 
closer  to  the  3-D  solution.  In  the  case  6  where  the  winding  angles  are  linearly  changing 
along  the  radius,  the  3-D  solutions  at  both  top  and  middle  surfaces  coincide  with  each 
other  and  can  be  best  predicted  by  the  PSS  solution. 

As  a  conclusion,  both  PSS  and  GPS  solutions  yield  the  better  solution  than  PSN  in 
all  cases.  Since  the  implementation  of  GPS  analysis  into  a  computer  code  takes  more 
efforts  than  the  PSS,  the  PSS  solution  is  the  best  in  prediction  the  stresses  and  strains  of 
the  flywheel  rotor,  regardless  of  the  hybridization  and  helical  winding,  and  the  thermal 
and  rotational  forces.  However,  care  should  be  taken  in  applying  the  2-D  analysis  to  the 
rotor  with  abrupt  changes  of  fiber  winding  angles. 
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a 


Figure  2.  One  eighth  of  composite  flywheel  rotor  for  3-D  finite  element  analysis. 


Table  1.  Specifications  of  the  cases 


Hybrid  Rotor 

Helical  Rotor 

Case 

1 

:  2 

3 

4 

5 

.  6 

ID  (mm) 

75 

75 

75 

75 

75 

75 

Thickness 

(mm) 

36.5/36.5 

36.5/36.5 

36.5/36.5 

36.5/36.5 

36.5/36.5 

73.0 

Height  (mm) 

35 

220 

35 

35 

220 

35 

Materials 

E- 

glass/Epoxy 

T300/5208 

E- 

glass/Epoxy 

T300/5208 

E- 

glass/Epoxy 

T300/5208 

T300/520 

8 

T300/520 

8 

T300/520 

8 

Angles  (  °) 

0/0 

0/0  | 

0/0 

45/0 

45/0 

45/.. 7  0 

AT  (°C) 

-100  ' 

-100 

• 

-100 

-100 

-100 

RPM 

• 

• 

10.000 

• 

• 

• 

Table  2.  Material  properties 


Material  data 

T300/5208 

E-glass/Epoxy 
(  v,  =  0.45  ) 

E-glass/Epoxy 
(  v ,  =  0.65) 

Ex  (GPa) 

181 

38.6 

56.96 

Ej=  E,  (GPa) 

10.3 

8.27 

8.27 

Vxy=Vxz  \ 

0.28 

0.26 

0.26 

Vyz  ? 

0.3 

0.3 

0.3 

Gxv  =  Gxz  (GPa) 

7.17 

4.14 

7.16 

Gyz(GPa) 

3.4 

3.4 

3.4 

ax  (pm/m)/K 

2.0  E-6 

8.6  E-6 

6.12  E-6 

ay=az  (pm/m)/K 

22.5  E-6 

2.1  E-6 

16.2  E-6 

p  (kg  /m3) 

1600 

1800 

2066 

15 


radius  [m] 

(a)  radial  stresses 


radius  [m] 

(b)  hoop  stresses 


0.0375  0.0475  0.0575  0.0675  0.0775  0.0875  0.0975  0.1075 
radius  [m] 

(d)  radial  strains 


radius  [m] 


(e)  hoop  strains 


radius  [m]  radius  [m] 

(c)  axial  stresses  (f)  axial  strains 


Figure  4.  Distribution  of  radial,  hoop,  and  axial  stresses  and  stains  of  case  2 
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5.  PROCESS  INDUCED  RESIDUAL  STRESSES 


The  flywheel  rotor  of  composite  materials  has  characteristics  of  distinctively  high 
energy  density,  long  life  and  lightweight.  An  advanced  composite  flywheel  has  thus 
attracted  the  considerable  attentions  of  many  investigators  for  various  energy  storage 
applications  such  as  electric  utilities,  hybrid  or  electric  vehicles  and  spacecraft  [1-5]. 
Recent  efforts  in  development  of  the  flywheel  have  been  devoted  to  the  experimental  and 
theoretical  investigation  on  increasing  the  performance,  while  insuring  the  safety  of  the 
flywheel,  along  with  technical  innovations  in  high  stiffness  and  strength  composite 
materials. 

In  the  design  of  such  high  performance  flywheel  rotor,  how  to  hybridize  the 
composite  materials,  i.e.,  determination  of  the  material  sequence  and  the  thickness,  is 
crucial  in  minimizing  the  radial  stresses  caused  by  the  high  rotating  speed  [6].  The 
influences  of  the  centrifugal  forces  on  the  mechanical  behavior  of  the  rotor  have  thus 
been  successfully  analyzed  in  design  of  a  thick  wound  composite  rotor  [7-10].  However, 
as  the  thickness  of  the  rotor  increases,  manufacturing  process  of  the  thick  section  polymer 
composite  rotors  can  generate  the  residual  stresses  which  add  to  those  generated  by  the 
centrifugal  forces  and  significantly  deteriorate  the  performance  of  the  flywheel  rotor 
[10,11]. 

Wet  winding,  followed  by  the  stages  of  heat  buildup,  curing,  and  cooling,  is  a 
typical  process  for  manufacturing  the  thick  composite  rotor.  During  heat  buildup,  the 
radial  stiffness  decrease  and  the  stresses  are  relaxed,  and  upon  curing  the  radial  stiffness 
and  transverse  strength  increase  and  simultaneous  chemical  shrinkage  occurs.  Cooling  is 
associated  with  simultaneous  thermal  shrinkage  of  both  the  rotor  and  the  mandrel  [12-16]. 
Especially,  the  cool-down  stage  in  the  curing  process  cause  the  tremendous  amount  of 
radial  tensile  stresses  primarily  due  to  the  anisotropic  thermal  expansion  of  individual 
plies.  During  the  entire  history  of  the  manufacturing  process,  the  mechanical  properties 
and  stress- strain  state  may  undergo  a  significant  change. 

Several  investigators  have  developed  theoretical  models  of  the  whole  filament 
winding  process  [13-17].  Among  several  attempts,  Lee  et  al.  modeled  the  processes 
applicable  to  the  moderately  thick  cylinders  [14-16]  by  establishing  submodels,  each  for 
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thermochemical  behavior,  fiber  motion,  void,  and  stress,  and  developed  a  cure  simulation 
computer  program  called  WINDTHICK  [16].  Even  though  the  models  are  well 
established,  we  should  supply  many  parameters  for  the  material  properties  to  simulate  the 
program.  Among  them  the  material  properties  of  the  uncured  composite  fiber  and  resin 
and  their  variation  along  with  the  temperature  and  the  degree  of  cure  are  rarely  handy. 
Therefore,  the  complexities  involved  in  these  models  makes  it  difficult  to  be  used  as  a 
design  tool,  especially  for  the  optimization  of  the  composite  cylinder. 

For  a  design  tool  for  the  thick  wound  cylinder,  a  relatively  simplified  model  was 
proposed  which  mainly  consider  the  thermal  expansion  and  detachment  of  the  mandrel 
and  the  cooldown  of  the  cylinder  [11].  In  the  work  it  is  found  that  main  factors 
influencing  the  stress  state  for  some  helically  wound  composite  cylinders  are  the 
temperature  change  and  the  chemical  shrinkage,  and  not  the  fiber  tension  nor  motion.  In 
that  model,  however,  the  detachment  between  the  mandrel  and  the  cylinder  should  be  still 
modeled,  and  the  material  properties  of  the  uncured  composites  need  to  be  given.  It  is 
observed  in  other  work  that  the  thickness  and  the  material  properties  of  the  mandrel  can 
affect  the  residual  stresses  in  the  composite  cylinder  [19].  However,  those  stresses  vary 
according  to  the  thickness  ratio  of  the  mandrel  and  the  composite  cylinder. 

It  has  been  observed  that  the  curing  process  for  moderately  thick  (e.g.  2”  thickness 
for  4”  ID)  graphite  wound  composite  cylinder  can  even  cause  a  radial  delamination  [11]. 
However,  the  experimental  measurement  of  the  residual  stresses  and  strains  of  the  thick 
wound  rotor  has  been  rarely  performed  in  comparison  to  the  relatively  thin  rotors  [18-20]. 
Moreover,  the  detection  of  the  delamination  present  in  the  thick  wound  cylinders  never 
been  reported. 

Compared  to  the  conventional  one-stage  cure,  the  residual  stress  advantages  can  be 
realized  by  the  multi-staged  curing  method,  where  layers  of  materials  are  allowed  to 
harden  for  a  time  period  before  more  materials  are  wound  onto  the  rotor.  An  in-situ 
curing  can  be  a  more  productive  way  than  the  staged  curing  method,  and  a  simplified 
analysis  of  residual  stresses', related  to  this  process  has  been  investigated  by  Bakis  et  al. 
[21].  Nevertheless,  the  one-stage  cure  is  always  a  more  productive  method  as  long  as  the 
residual  stresses  are  properly  managed. 

In  this  study,  in  order  to  estimate  the  stresses  and  strains  accumulated  inside  the 
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rotor  during  the  manufacturing  process  of  the  thick  wound  cylinder,  a  simplified 
theoretical  analysis  using  the  material  properties  only  at  room  temperature  and  the 
equivalent  thermal  loading  is  presented.  The  validity  of  the  method  is  examined  by 
comparison  with  the  results  from  the  cure  simulation  program  WINDTHICK.  The 
experimental  test  was  also  performed  where  the  residual  stresses  were  released  by 
splitting  the  rotor  in  the  radial  direction,  while  measuring  the  strains  [22,23],  Most 
importantly,  the  residual  stresses  and  strains  that  are  high  enough  to  cause  the  radial 
delamination  have  been  studied  and  experimentally  measured. 


5.1  Stress  and  Strength  Analysis  of  a  Hoop  Wound  Cylinder 


A  stress  analysis  of  a  hoop  wound  composite  cylinder  is  performed  with  a 

* 

consideration  of  the  process  induced  residual  strains  that  mainly  result  from  the 
temperature  drop  and  the  chemical  shrinkage.  A  linear  plane  stress-strain  relationship 
considering  the  thermal  strains  and  chemical  shrinkage  in  a  cylindrical  coordinate  system 
can  be  written  as 
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where  Gr ,  Ge,  £r  and  £e  are  the  radial  and  hoop  stresses  and  strains,  respectively. 
The  coefficients  of  thermal  expansion  (CTE)  are  denoted  by  ae  and  ar .  A  typical 


temperature  change  A Tc  in  the  curing  process  is  considered  as  the  difference  between 
the  maximum  cure  temperature  and  the  room  temperature.  In  general,  the  chemical 
shrinkage  occurring  in  the  process  can  be  hardly  measured.  However,  since  both  thermal 
strains  and  chemical  shrinkages  of  the  composite  materials  mainly  result  from  the 
deformation  of  the  matrix,  each  magnitude  of  the  chemical  shrinkage  in  the  hoop  and 

radial  directions  is  assumed  to  have  the  same  ratio  as  the  thermal  strains.  Hence,  Eq.  (34) 
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can  be  written  as  \ 


Qn  Qn 
Ql\  Ql2. 


fo  \ 


\ErJ 


( a,  ^ 


AT 


(35) 


where  an  effective  temperature  change  AT  contains  both  ATc  and  a  temperature  drop 
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A Tsh  equivalent  to  the  chemical  shrinkage: 


AT  -  AT  +  AT. 


(36) 


To  derive  the  stresses  and  strains  due  to  the  effective  temperature  change  A T ,  Eq.(35)  is 
rewritten  as 
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where  p  =  Qa  . 

The  stress  distribution  of  a  hoop  wound  composite  cylinder  is  governed  by  the 
radial  equilibrium  equation,  which  is  written  in  cylindrical  coordinates  as  [24-26] 


gfr  |  <rr-0e  -0 
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The  circumferential  and  radial  strains  are  expressed  by  the  r-direction  displacement  ur : 


ur  dur 

*e= —>£r=^r 
r  dr 


(39) 


Substituting  Eqs.  (37)  and  (39)  into  Eq.  (38)  yields  the  radial  and  hoop  stresses  and 
strains: 


ar  =  +  C2(p17~K~1  +  (pT1Af 

(40a) 

£r  =  xCj^F*"1  +  KC2(p^~K~'  +  <pnAT 

(40b) 

<7e  =  xC,r ,'1  -  kC27~*~]  +  (pT2AT 

(41a) 

£e  —  C^r  —  C2(p2r  +(pnAf 

(41b) 

where  the  radial  coordinate  is  normalized  by  the  inner  radius  r„  denoted  by  r  .  The 
symbols  /rand  (pt  are  defined  in  terms  of  the  material  properties: 


J^-,  (p{  = - - - ,  (p2~ - - - , 

W*  1  *222+02  *222-02 

*Pt\  =  n-_nr  9  anC^  tpn  =  (2l2  +  022  )  fin  —  fir 
^22 


(42a). 

(42b) 


Unknown  constants  C\  and  Ci  are  determined  from  the  boundary  conditions,  i.e.  ar= 0  at 
r=Vi  and  r=r0  : 
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where  r0  is  the  dimensionless  outer  radius  normalized  by  the  inner  radius. 

The  radial  delamination  can  occurs  within  the  cylinder  due  to  the  residual  stresses. 
In  the  absence  of  the  shear  stress,  two  simple  delamination  criterions  can  be  used,  i.e.,  a 
maximum  stress  criterion  and  a  maximum  strain  criterion  [25]: 
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Delamination  occurs  when  the  delamination  index  Da  or  De  reaches  the  value  of  one. 
Once  a  delamination  occurs,  it  is  assumed  that  the  cylinder  is  split  into  two  cylinders. 

The  radial  location  where  the  maximum  delamination  index  occurs  is  derived  from  the 


condition  that  =  0  in  Eq.  (44),  each  yielding 
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It  is  noted  that  the  radial  location  a ,  where  the  maximum  index  Dc ,  or  the 


maximum  radial  stress  occurs,  is  a  function  of  only  T0  and  k,  independent  of  the  CTE 
a#  and  ar .  And  7^  £  is  always  located  farther  than  7^  a  from  the  center  of  the 

cylinder.  Substitution  of  Eq.(45a)  into  Eq.(40a)  then  yields  the  maximum  radial  stress; 
Eq.(45b)  into  Eq.(40b)  yields  the  maximum  radial  strain.  Using  these  formulae  together 
with  either  Eq.(44a)  or  (44b),  the  maximum  temperature  drops  which  cause  delamination 
are  obtained,  each  denoted  by  ATd  c  and  ATd  E  : 


A  Td_a  = 


\ _ (C-r'V  _ 


(46a) 


26 


{r^-zT^llA 


<pt2  a 


■Ofi+fn'ft 


(46b) 


It  is  noted  that  both  ATd  a  and  A Td  £  are  functions  of  the  dimensionless  radius 
T0  =  r0  /  T; .  In  other  words,  the  cylinders  of  the  geometric  similarity  experience  a 
delamination  at  the  same  temperature  drop.  It  can  be  thus  deduced  that  the  larger  cylinder 
is  more  advantageous  in  decreasing  the  thermal  stresses,  especially  in  the  flywheel  rotor. 


5.2  Cure  Simulation  by  WINDTHICK 

The  cure-simulation  program  WINDTHICK  is  used  to  determine  the  major  factors 
affecting  the  residual  stresses  within  the  thick  wound  composite  cylinder  among  the 
process  parameters  or  models,  which  include  (1)  fiber  winding  tension  and  motion,  (2) 
thermochemical  evolution,  (3)  interaction  with  the  mandrel,  and  (4)  change  of  composite 
material  properties  with  respect  to  the  temperature  and  the  degree  of  cure.  The  simple 
stress  analysis  presented  in  the  previous  chapter  is  then  performed  using  only  the  material 
properties  of  the  cured  composites  rather  than  considering  in  detail  the  variation  of  the 
material  properties  along  with  the  degree  of  cure.  It  is  then  shown  that  the  total  residual 
stresses  can  be  well  represented  by  the  effective  temperature  drop. 

To  validate  the  approach,  residual  stresses  are  calculated  for  thick  cylinders  in  two 
cases  which  was  considered  by  Lee  et  at.  [16].  Fiberite  T300/976  is  used,  and  the  material 
properties  and  cure  cycles  are  summarized  in  [16].  Composite  cylinders  are  wound  on  a  4 
inch  diameter,  0.25  inch  thick  aluminum  mandrel.  Two  stacking  sequences  are  considered, 
i.e.,  [-55/+55]so  and  [(±45)I5,9040,(±45)15] .  Each  cylinder  consists  of  100  layers,  and  the 
corresponding  wall  thickness  is  2  inches.  Room  temperature  and  maximum  cure 
temperature  are  taken  respectively  as  65  °F  and  350  °F  [16].  The  calculation  is  performed 

in  two  ways;  one  with  the  full  model  considered  by  Lee  et  al.  [16]  and  the  other  with  the 

\ 

simple  model  in  which  the  cured  material  properties  are  used  together  with  the 
temperature  drop  and  the  chemical  shrinkage.  In  the  simple  model,  the  temperature  drop 
for  both  cases  is  given  as  the  difference  between  the  room  and  maximum  cure 
temperatures,  i.e.  Tc=  65  °F  -350  °F  =  -285  °F,  and  the  chemical  shrinkage  as  £sh  =0.5%. 
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The  calculated  radial  and  hoop  stresses  are  in  good  agreement  with  the  results  by  the 
program  WINDTHICK  within  3%  of  error. 

As  a  second  case,  the  effects  of  the  rotor  thickness  on  the  residual  stresses  are 
investigated.  Three  cylinders  A,  B  and  C  of  different  outer  radii,  i.e.  77,  1 10.5,  and  152.5 
mm  are  considered  with  the  same  37.5  mm  inner  diameters,  and  they  are  all  hoop  wound. 
The  cure  cycle  is  chosen  to  fully  cure  the  cylinders.  All  other  parameters  such  as  the 
mechanical  properties  of  the  mandrel  and  the  composites  are  the  same  as  the  first  case. 
The  radial  and  hoop  stresses  are  calculated  and  shown  in  Figure  9.  It  is  once  again  shown 
that  the  results  by  the  simple  stress  analysis  are  very  close  to  those  by  the  full  model.  It  is 
noticed  that  the  same  temperature  drop  of  Tc=  -285  °F  and  chemical  shrinkage  of  0.5  % 
for  all  three  cylinders  in  the  simple  model  can  represent  the  residual  stresses  from  the  full 
model. 

It  is  concluded  in  this  section  that  the  temperature  drop  and  the  chemical  shrinkage 
are  the  major  sources  for  the  residual  stresses  within  the  thick  wound  composite  cylinder, 
not  the  interaction  with  the  mandrel  nor  the  fiber  tension  and  motion.  In  other  words,  the 
distribution  of  the  process-induced  residual  stresses  can  be  obtained,  regardless  of  the 
rotor  thickness,  from  the  cured  composite  cylinder  experiencing  a  uniform  temperature 
drop  and  chemical  shrinkage  throughout  the  rotor.  The  results  by  the  simple  stress 
analysis  can  be  also  obtainable  from  the  WINDTHICK  program  with  the  exclusion  of  the 
thermal  expansion  of  the  mandrel,  the  fiber  tension  and  motion  and  the  dependency  of  the 
material  properties  upon  the  degree  of  cure. 


\ 
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5.3  Split  Ring  Test  Method 


The  residual  strains  are  measured  using  a  split  ring  method.  The  strain  gages  are 
glued  to  the  side  surfaces  of  the  cured  cylinder.  The  cylinder  is  then  split  by  cutting  along 
the  radial  direction,  and  the  residual  stresses  are  released.  The  stress-released  strains  are 
measured. 

Assuming  that  the  stresses  and  strains  are  zero  through  out  the  rotor  just  before 
cooling  down  from  the  cure  temperature,  the  residual  strains  and  stresses  have  been 
accumulated  during  the  effective  temperature  drop  A  f .  Theses  residual  stress  and  strain 
state  can  be  viewed  as  a  sum  of  two  consecutive  states;  stress-free  temperature  drop 
followed  by  application  of  the  mechanical  stresses  to  the  radial  surfaces.  The  typical 
stress  and  strain  distribution  and  a  superposition  procedure  are  illustrated  in  Figure  10. 
The  stains  of  the  stress-free  specimen  due  to  the  pure  temperature  drop  are  simply  written 
as 

eT=aAf  (47) 

On  the  other  hand,  the  relationship  of  the  mechanical  stresses  and  strains  generated  by 
applying  only  the  mechanical  loading  to  the  surface  can  be  written  as 

gm  =  Q$m  (48) 

The  residual  stresses  and  strains  of  the  full-circled  rotor,  each  denoted  by  (Jr  and  £r,  are 
then  obtained  by  superposing  stresses  and  strains  in  Eqs.  (47)  and  (48)  and  : 

cR  =  cr  M ,  er  —  eM  +  a  AT  (49) 

Since  the  split  of  the  specimen  causes  the  stress-free  boundary  conditions  at  the  cutting 
surface,  cutting  the  specimen  in  the  radial  direction  is  mechanically  equivalent  to 
applying  the  negative  of  the  circumferential  components  of  the  residual  stresses  to  the 
radial  surface.  In  other  words,  the  residual  strains  eM  are  released  and  measured  on 

mechanicaL  split  of  the  specimen. 

\ 

The  simple  but  effective  superposition  justifies  the  usefulness  of  the  split  method 
for  measuring  the  residual  strains  and  stresses  accumulated  in  the  thick-walled  rotor  made 
of  composite  materials.  Notice  that,  in  the  case  of  the  hybrid  composite  rotor,  the  stresses 
and  strains  of  the  free-edged  specimen  caused  by  the  temperature  drop  should  be 
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calculated  using  a  finite  element  method  since  this  case  is  no  longer  axisymmetric.  For 
completeness,  each  step  is  verified  using  three-dimensional  finite  element  methods.  The 
total  residual  stresses  and  strains  can  be  always  calculated  by  the  stress  analysis  described 
in  the  previous  section. 


\ 


Figure  10.  The  residual  stresses  and  strains  developed  inside  the  circular  rotor  are 
superposition  of  those  by  the  stress-free  temperature  drop  plus  those  by  the  mechanical 
loadings  applied  to  the  cut  surface. 
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5.4  Measurement  of  Residual  Strains 


The  composite  materials  are  ideal  for  the  flywheel  rotor  since  the  stresses  generated 
by  the  rotational  centrifugal  forces  can  be  effectively  supported  by  the  anisotropic 
materials.  On  the  other  hand,  however,  the  anisotropic  coefficients  of  thermal  expansion 
differing  in  the  radial  and  circumferential  directions  inevitably  generate  the  residual 
strains  proportionally  to  the  temperature  drop.  In  addition,  the  distribution  of  these 
residual  strains  is  affected  by  the  dimensional  size  and  the  stiffness  anisotropy  These  can 
be  well  observed  in  the  previous  section.  Most  importantly,  understanding  size  effects  on 
the  residual  strains  is  critical  in  design  of  the  flywheel  rotor  made  of  composite  materials. 

To  validate  the  effective  temperature  change  presented  in  the  previous,  section,  the 
residual  strains  were  experimentally  measured  for  three  different  outer  radii  of  77,  110.5, 
152.5  mm  with  the  same  inner  radii  of  37.5  mm  and  the  length  of  25.4  mm;  each 
specimen  is  denoted  respectively  by  specimen  A,  B,  and  C.  The  aluminum  mandrel  of 
outer  and  inner  radii  of  32.5  mm  and  27.5  mm  respectively  were  first  machined  with  two 
side  circular  plates  of  the  outer  radius  of  150  mm  and  the  thickness  of  10  mm  attached  for 
guiding  the  winding  fiber  to  produce  the  rotor  of  length  of  200  mm,  as  shown  Figure  1 1 . 
Resin  system  was  prepared  by  mixing  epoxy,  catalysis  and  accelerator  with  the  mass 
fraction  ratio  of  60:35:5.  During  the  filament  winding  process,  graphite  Fibers  of  T300- 
6K,  followed  by  dipping  in  the  resin  bath,  were  wound  with  the  fiber  tension  of  10  N  onto 
the  mandrel  rotating  with  the  revolution  speed  of  50  rpm.  The  nose  which  fed  the  fibers 
axially  traveled  at  the  speed  of  lm/min.  Filament  winding  process  was  then  followed  by 
curing  in  the  autoclave  with  the  cure  cycle  as  shown  in  Figure  12.  In  order  to  ensure  the 
axisymmetric  material  properties,  the  one  end  of  the  horizontal  mandrel  shaft  was 
connected  to  the  motor  and  revolved  with  the  rotational  speed  of  15  rpm  during  the  whole 
curing  process.  Once  cure  was  completed,  the  rotor  was  spliced  and  machined,  yielding 
the  length  of  25.4  mm.  i 

In  order  to  assure  the  quality  of  the  cylinder,  the  fiber  volume  fractions  are 
measured  along  the  radial  direction  as  shown  in  Figure  13a.  The  distribution  of  the  fiber 
volume  fractions  along  the  radial  direction  are  shown  in  Figure  13b,  and  the  average 
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fraction  is  68%.  Photomicrographs  of  the  cross-sections  are  shown  in  Figure  13c. 

Two  sets  of  strain  gages  in  the  radial  direction,  each  for  radial  and  circumferential 
strains,  were  then  glued  on  the  one  side  surface  of  the  rotor  as  shown  in  Figure  14.  Each 
set  consisted  of  from  3  to  15  stain  gages  depending  on  the  size  of  the  rotors.  The  residual 
stresses  were  released  by  the  split  ring  method,  that  is,  cutting  D-section  out  of  the  rotor. 
The  readouts  before  cutting  the  D-section  were  set  to  zero,,  and  then  recorded  again  after 
cutting.  Note  that  the  mechanical  residual  strains  are  opposite  in  sign  of  those  recorded 
upon  cutting.  These  values  from  the  strain  gages  for  each  size  of  rotor  are  shown  in 
Figure  15:  the  measured  radial  and  circumferential  strains  of  the  specimens  A,  B,  and  C 
are  represented  by  solid  dots  as  shown  in  (a)  through  (f),  respectively.  Whereas  the  radial 
and  circumferential  strains  of  the  specimen  A  are  continuously  distributed  along  the 
radius,  the  abrupt  changes  are  characteristic  of  the  strains  of  the  specimens  B  and  C.  The 
abrupt  change  of  the  strains  along  the  radius  was  not  clearly  understood  until  the  stress 
analysis  for  the  residual  strains  was  performed. 

The  effective  temperature  change  is  found  to  be  140  °C  to  best  fit  the  experimental 
data  in  the  specimen  A  as  shown  in  Figure  15.  The  material  properties  used  in  the 
calculation  are  from  those  for  T300/Epoxy.  It  is  found  that  the  temperature  drops  less 
than  140  eC  for  the  specimens  B  and  C  cause  the  delamination  in  the  radial  direction. 
This  is  manifested  in  Figure  16:  the  radial  and  circumferential  stresses  of  the  specimen  A, 
B,  and  C  with  the  temperature  drop  of  140  eC  are  calculated  and  respectively  shown  in 
solid  lines  in  Figure  16(a)~(f).  It  can  be  seen  from  Figure  16(c)~(f)  that  the  strength  ratios 
for  the  specimens  B  and  C  are  beyond  the  value  of  1  with  the  temperature  drop  of  140°C , 
which  means  that  the  delamination  have  occurred  before  the  temperature  drop  of  140°C . 
The  temperature  drops  of  93  and  82  °C  are  calculated  which  initiate  the  delamination  of 
the  specimens  B  and  C,  respectively.  The  radial  and  circumferential  stresses  at  these 
temperature  drops  are  shown  in  dotted  lines  in  Figure  16(c)~(f).  After  these  delaminating 
temperature  drops,  the  specimen  are  assumed  to  be  split  into  two  rings,  and  a  further 
temperature  is  dropped,  thus  totally  140  CC.  These  consecutive  residual  strains  of  the 
specimen  B  and  C  are  shown  in  solid  lines  in  Figure  15(c)~(f).  It  is  noticed  that  the 
damaged  area  experience  the  high  radial  strains  and  stepped  aside  the  circumferential 
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strains.  Even  though  the  locations  of  the  delamination  predicted  for  the  specimens  B  and 
C  do  not  exactly  match  with  those  measured  by  the  split-ring  test,  overall  distribution  of 
the  calculated  residual  strains  are  in  good  agreement  with  the  measured  ones. 

As  a  conclusion  in  this  chapter,  the  residual  stresses  generated  during  the  elevated 
temperature  cure  of  thick  section  polymer  composite  components  can  lead  to 
delamination.  The  purpose  of  this  research  was  to  gain  a  better  understanding  of  residual 
strains  in  a  simple  way,  and  develop  methods  for  measuring  and  predicting  them.  These 
could  then  be  used  to  investigate  ways  to  manage  residual  stresses  in  order  to  reduce  the 
residual  stresses  and  stains.  Including  the  contact  between  the  mandrel  and  rotor  does 
change  the  results,  but  low  enough  to  be  neglected. 


Figure  11.  Mandrel  with  two  guide  plates  for  filament  winding 


\ 

t 


34 


Figure  13.  Measurement  of  the  fiber  volume  fraction  along  the  radial  direction 
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(e)  Specimen  C  before  split  (f )  Specimen  C  after  split 


Figure  14.  Three  specimens  A,  B  and  C  before  and  after  split  in  the  split-ring  test; 
the  stresses-released  strains  are  measured. 


36 


5.5  Summary 

The  process-induced  residual  strains  can  lead  to  delamination  within  a  thick  wound 
composite  cylinder.  Especially  in  the  case  of  the  flywheel  rotor,  the  residual  radial 
stresses  add  to  those  generated  by  the  centrifugal  forces,  which  significantly  deteriorate 
the  performance  of  the  rotor  since  the  radial  tensile  stresses  are  mostly  sustained  by  the 
relatively  weak  matrix.  The  process-induced  residual  strains  should  be  thus  properly 
considered  in  design  of  the  thick  wound  composite  cylinders.  However,  a  detailed 
universal  rheological  model  including  cure  shrinkage  as  well  as  thermal  strains,  the 
process  of  gelation,  cure  kinetics,  the  exotherm  and  interaction  between  the  rotor  and  the 
mandrel  would  be  very  difficult  to  set  up  for  design  of  the  thick  wound  composite  rotor. 

The  program  WINDTHICK  has  been  used  to  generate  the  residual  strains  and  it  has 

$ 

been  shown  that  the  main  factors  influencing  the  residual  stresses  in  the  thick  wound 
cylinder  are  mainly  due  to  the  temperature  change  and  the  chemical  shrinkage.  A  simple 
stress  analysis  with  an  effective  temperature  change  is  then  proposed  which  efficiently 
predict  the  process-induced  residual  stresses.  In  this  way,  the  stress  analysis  can  be 
effectively  linked  to  an  optimal  design  of  the  thick  wound  cylinders  with  a  consideration 
of  the  process  induced  residual  strains.  Split  ring  tests  have  been  performed  to  measure 
the  residual  strains  within  the  thick  wound  cylinders  of  different  sizes.  The  measured  data 
are  in  good  agreement  with  the  results  by  the  present  analysis. 


6.  OPTIMIZATION  COMPUTER  CODE  -  FLYOPT 


A  proper  combination  of  composite  materials  for  the  rotor  could  reduce  the  stress 
distribution  in  the  rotor  and  eventually  increase  TSE.  The  objective  of  this  study  is  to 
maximize  TSE  of  a  hybrid  flywheel  rotor  without  any  material  failure.  Strength  ratio  R  is 
used  for  failure  analysis,  where  a  value  of  R  greater  than  1  indicates  material  failure  [8]. 
To  ensure  the  rotor  against  any  material  failure,  in  this  approach,  strength  ratio  R  and  the 
damage  indices  Dc  and  De  calculated  at  any  point  of  the  rotor  must  have  values  less  than 

1 .  The  inner  radius  r\  and  height  h  of  the  rotor  are  given  and  the  thickness  tm  of  each  rim 
(w=l,2,3,...,Af)  is  used  as  the  design  variables.  The  optimization  is  then  formulated  as 
follows: 


Find  tm  (w  =  1,2, 3,-, A/) 

1  9 

Maximize  E  ~~Ico2  (50) 

Subjectto  MaxR,DofD£  <1 

where  I  is  a  mass  moment  of  inertia  of  the  hybrid  rotor,  and  CO  is  a  rotating  speed.  The 
total  stored  energy  (TSE)  E  can  be  thus  written  in  terms  of  the  inner  radius  rm  and  outer 
radius  rm+j  of  each  m- th  rim: 


E  =  !7  ha>2fJ{rL~r:)pm  (51) 

4 

where  the  subscript  m  indicates  each  rim  and  the  other  variables  include  the  length  of  the 

rotor  h ,  and  the  density  of  the  m- th  rim  pm . 

In  order  to  solve  the  nonlinear  optimization  problem  formulated  in  Eq.  (50),  the 
\ 

modified  method  of  feasible  directions  for  constrained  minimization  is  used.  The  main 
task  in  the  optimization  method  is  to  find  a  usable-feasible  search  direction  and  the  one¬ 
dimensional  search.  In  the  every  iteration  of  optimization,  stress  and  strength  analysis  are 
performed  and  the  maximum  strength  ratio  in  each  ring  is  evaluated  as  described  below. 
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An  optimal  design  of  a  hybrid  composite  rotor,  by  maximizing  TSE  without  any 
material  failure  as  formulated  in  Eq.  (50),  is  performed  using  a  combination  of  composite 
materials  with  different  fiber  types:  (A)  Glass/Epoxy,  (B)  T300/2500,  (C)  T800H/2500. 
Optimizations  are  performed  combining  up  to  three  different  composite  rims  and  can  thus 
be  grouped  into  3  categories.  In  the  first  category,  only  one  material  is  used,  and  in  the 
second  category,  two  materials,  and  so  on.  In  each  category  only  several  cases  of  material 
selections  and  sequences  are  considered  because  it  is  efficient  way  of  increasing  TSE  to 
arrange  softer  materials  in  the  inner  side  and  stiffer  materials  in  the  outer  side  of  the  rotor. 
The  thickness  of  each  composite  rim  is  considered  as  a  design  variable  and  equally 
divided  into  element  rings,  resulting  in  a  total  of  30  element  rings.  The  initial  thickness 
for  each  composite  rim  is  20  mm,  and  the  inner  radius  of  the  rotor  is  50mm.  The  value  for 
the  rotating  speed  is  taken  as  60,000  RPM. 

All  the  calculated  optimum  TSE  are  shown  in  Figure  17  where  each  bar  indicates 
the  location  of  each  rim  and  thus  the  length  of  the  bar  means  the  thickness  of  each  rim.  As 
noted  in  the  stress  distribution,  the  effects  of  the  curing  temperature  on  the  optimal  design 
are  significant.  For  instance,  the  results  of  the  optimal  designs  using  the  material 
sequence  A-B-C  show  that  TSE  reduces  by  about  30%  with  AT  =  -110°C  and  36%  with 
A  f~  -125  °C.  In  other  words,  neglecting  the  residual  stresses  may  leads  to  causing  the 
higher  stresses  than  expected  in  the  process  of  operating  the  flywheel  rotor.  The 
distributions  of  the  strength  ratios  for  A f  =  0°C,  AT-  -110°C  and  A  f-  -125°C  are 
shown  in  Figure  6.  It  is  noted  that  the  maximum  TSE  is  obtained  when  the  strength  ratio 
R=1  is  reached  in  every  composite  rim.  Thus,  strength  is  not  wasted  in  each  rim  which 
otherwise  experience  proportionately  smaller  stress  levels.  This  means  that  maximizing 
TSE  is  equivalent  to  incorporating  the  highest  possible  stresses  that  satisfy  the 
equilibrium  and  compatibility  conditions.  At  the  optimum  state,  each  thickness  of  the  rim 
is  thus  determined  so  that  the  strength  ratio  R=  1  is  reached  in  every  composite  rim.  It  can 
be  obviously  noted  that  the  number  and  sequence  of  materials  are  very  important  in 
maximizing  TSE  or  minimizing  the  stresses  inside  the  materials.  The  stress  distributions 
are  effected  by  the  sequence  of  stiffness,  mostly  hoop  directional  Young’s  modulus. 
Softer  materials  inside  of  the  rotor,  which  are  easier  to  expand  in  the  radial  direction,  and 
stiffer  materials  outside  of  the  rotor,  which  experience  less  radial  displacement, 
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eventually  reduce  the  radial  stresses. 

A  computer  program  FLYOPT  is  developed  using  Microsoft  Excel  user-interface. 
The  input  data  in  blue  and  the  output  data  in  red  for  the  stress  analysis  and  the  user- 
interface  for  the  optimization  are  respectively  shown  in  Figures  18  and  19.  The  design 
variables  in  the  optimization  module  can  be  linked  to  any  variables  in  the  stress  analysis 
module. 


0  25  50  75  100  125150  175  20®  25  50  75  100  125  150  175  20(0  25  50  75  100125150175200 

radius,  mm  radius,  mm  radius,  mm 

(a)  Af=0°C  (b)Af=-110°C  (c)Af=-1250C 


Figure  17.  Comparison  of  optimal  thickness  of  composite  rims  and  the  corresponding  TSE 
(Wh)  for  various  selections  of  composite  material  (<2r=60,000RPM  ). 
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INPUTS  in  Blue 


w(RPM) 

65392 

Max  w  (RPM) 

32696 
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RJ: 

0.0095 

R__o : 

0.0762 

R_m: 

0.0095 

Total_Thickness 

0.0667 

Nel  : 

150 

plane-option;  0  for  PSr 

2 

Solution  Type 
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Hybrid  Rings 

Mandrel 

Group_J 

Group_2 

Group_3 

Group_4 

Material  Group  Numb 

0 

1 

2 

3 

4 

11 

11 

11 

11 

11 

|Matcrial_Thickncss 

0 

0.00635 

0.009525 

0.009525 

E5E3ES3MB 

Scotch] SI);  vf=0.65 

otch|SI]:  vf=0.65 

tchJSIj:  vf-0.65 

>tchfSI]:  vf=0.65 

:>tch|SIJ;  vN0.65 

|Number_of_Rings 

0 

20 

20 

20 

20 

0 

0 

0 

0 

0 

Material_temp_Currei 

25 

25 

25 

25 

25 

Material_dT 

25 

25 

25 

25 

25 

Matcrial_chemical_sh 

0 

0 

0 

0 

0 

Material_Group 

0 

20 

40 

60 

80 

RESULTS  in  Red 


Ip  (Kg* m2) 

QE3EC3IIBH 

0.94 

Volume  (liter) 

0.456093736 

Ejnax  (Wh) 

18 

5 

E_usable 

14 

SED_max  (Wh/Kg) 

19.20 

SEDjmin(Wh/K[ 

4.80 

SED_usable(Why 

14.4028545 

SVED_max  (Wh/liter 

w 

VO 

©■> 

00 

SVED_max(Wh/l 

9.92 

29.75629739 

Tip  Speed  (m/sec) 

|  521.8049681 | 

BSEESEnS®. 

65.22562101 

|Cost material  (Sj 

2.768032882 

Failure  Analysis 

Sig_theta/X,MAX 

Sig_theta/X.MIN 

Sig_z/Y,  Max 

Sig_r/Y,  Max 

Tsai-Wu 

Winding 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00263 

-0.00567 

-0.00570 
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Figure  18.  Input  data  and  output  results  in  the  stress  module  of  the  program  FLYOPT  for 
the  flywheel  rotor  subjected  to  both  the  thermal  changes  and  the  rotational  speed. 
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Figure  19.  The  design  variables,  the  constraints,  and  the  objective  functions  specified  in 
the  optimization  module  of  the  program  FLYOPT. 
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7 .  CONCLUSIONS 


The  goal  of  the  present  work  is  to  develop  an  optimization  tool  for  a  design  of  the 
flywheel  rotor  with  a  consideration  of  the  process-induced  residual  strains  accumulated 
during  the  winding  and  curing  process.  A  stress  analysis  of  a  helically  wound  hybrid  rotor 
is  developed  in  which  an  equivalent  temperature  change  is  used  to  consider  the  residual 
stresses  and  strains.  The  validity  of  the  stress  analysis  has  been  proved  using  a  three- 
dimensional  finite  element  program.  The  effects  of  the  material  hybridization  and  the 
helical  winding  angles  oa  the  thermal  stresses  are  also  investigated.  A  WINDTHICK 
program  is  also  used  to  validate  an  effective  temperature  change  considered  in  the 
analysis.  The  residual  stresses  have  been  experimentally  measured  using  a  split-ring  test. 
As  a  product  of  this  work,  a  computer  program  for  the  optimum  design  of  the  hybrid 
composite  flywheel  rotor  are  developed  which  covers  the  various  size  and  material 
properties  of  the  rotor. 
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